Comments on the Redshift Distribution of 44,200 SDSS Quasars: 
Evidence for Predicted Preferred Redshifts? 
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ABSTRACT 

A Sloan Digital Sky Survey (SDSS) source sample containing 44,200 quasar redshifts is ex- 
amined. Although arguments have been put forth to explain some of the structure observed in 
the redshift distribution, it is argued here that this structure may just as easily be explained by 
the presence of previously predicted preferred redshifts. 

Subject headings: galaxies: active - galaxies: distances and redshifts - galaxies: quasars: general 



1. Introduction. 

It was realized early in the SDSS survey 
(Richards et al. 2002) that quasar redshifts near 
z = 2.7 would be under-sampled simply because 
at that redshift the colors of stars resemble those 
of quasars and it would be difficult to keep the 
quasar identification efficiency high. Based on the 
knowledge gained during the commissioning pe- 
riod, and the results from the first 2555 quasars, 
the target selection algorithm was refined to the 
point where it was felt by 2000 November that it 
was appropriate for use in obtaining a statistical 
sample of quasars (Richards et al. 2002). How- 
ever, after additional observations were carried 
out, and many more sources with redshifts above 
z = 3 were obtained, it became apparent that there 
was a second dip in the redshift distribution at z = 
3.5. This was again deemed to be a problem with 
the source selection algorithm, although it is not 
clear how this was determined, since a real source 
deficiency at this redshift would presumably give 
the same result. It may, however, have been de- 
cided from the evidence for a trough visible in 
the i-redshift plot (Schneider et al. 2003, Fig 3) 
which, when redshifts are assumed to be purely 
cosmological, can only be explained by a selec- 
tion effect. On the other hand, it is exactly what 
would be expected if the redshifts are intrinsic. 
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However, further changes to the target selection 
algorithm were made and it was finalized a second 
time on 2001 August 24. Spectra taken since that 
date were deemed to constitute the basis for the 
statistical sample of SDSS quasars (Richards et 
al. 2002). The first release of SDSS quasar red- 
shift data (DR1) contained the 16,713 redshifts 
measured prior to 2001 October (Schneider et al. 
2003). Although the target selection algorithm 
was being refined throughout most of the period 
the DR1 observations were made, the sample con- 
sidered here contains an additional ~ 27, 500 red- 
shifts that were obtained since 2001 October, after 
the algorithm was finalized. 

The primary SDSS quasar science goals were to 
study (1), the evolution of the quasar luminosity 
function and (2), the spatial clustering of quasars, 
as a function of redshift. However, a reasonably 
complete statistical sample of quasar redshifts can 
also be used to examine the long-standing contro- 
versy surrounding the very nature of quasar red- 
shifts. Peaks in the distribution of quasar redshifts 
have long been claimed as possible evidence for 
the existence of preferred redshifts (Karlsson 1971, 
1977; Burbidge and Napier 2001; Bell 2002a,c). 
However, only source samples with a few hundred 
redshifts have been available to test this claim. 
Now that surveys containing tens of thousands of 
quasar redshifts are available, if selection effects 
can be properly accounted for, although this will 
not be easy, any peaks still visible in the redshift 
distribution could be a much more reliable indica- 
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tor for the existence of preferred redshifts. Alter- 
natively, if predicted redshift peaks do not appear, 
it can mean either that intrinsic redshifts do not 
exist, or that they have been smeared out by a 
significant Dopplcr component. 

There are currently only two models that have 
been proposed that clearly define the preferred 
redshift values expected. It is of interest now to 
see if one of these models might equally well ex- 
plain the observed redshift distribution of 44,200 
SDSS quasars. 

2. Preferred Redshift Models 

Over 30 years ago it was pointed out (Karlsson 
1971, 1977) that the relation 

Alog(l+z) = 0.089 - — (1) 

appears to predict redshifts that are preferred 
in the emission line redshift distribution of the 574 
strong quasars detected in early surveys. Most of 
these redshifts were between z = and z = 2. Since 
that time much work has been carried out in an 
attempt to find further evidence to support this 
claim (Burbidge and Napier 2001). A significant 
amount of evidence has also been accumulated 
supporting the argument that high-redshift QSOs 
are ejected from nearby active galaxies, and that 
the redshifts of these objects are largely intrinsic 
(Arp 1998, 1999; Burbidge 1995, 1997a,b; Chu ct 
al. 1998; Pietsch et al. 1994; Lopez-Corredoira and 
Guticrcz 2002). Below z = 4.8, eqn 1 predicts pre- 
ferred redshifts of z = 0.061, 0.3, 0.6, 0.96, 1.41, 
1.96, 2.63, 3.45, and 4.47 (Burbidge and Napier 
2001). 

In addition to the model defined by eqn 1, a 
second intrinsic redshift model has been reported 
(Bell 2002a, c) that predicts intrinsic redshifts pe- 
riodic in z. This model resulted primarily from the 
study of the 15 high-rcdshift QSOs located around 
the Seyfert galaxy NGC 1068 (Bell 2002a,b,c; Bell 
and Comeau 2003a) which also assumed that the 
QSOs had been ejected from the Seyfert galaxy. 
In this model the high redshifts of the QSOs must 
contain large intrinsic components if they are at- 
tached to NGC 1068 which has a cosmological red- 
shift of z = 0.0038. 

From the combined results of the study of QSOs 
near NGC 1068, and the results obtained in an ear- 
lier one (Burbidge and Hewitt 1990), it was found 
that the intrinsic redshift components in quasars 



(zjq) can be defined by the relation 
z lQ = 0.62[N - M N ] - - (2) 

where N is an integer, and Mjy is a function of 
a second quantum number n, (Bell 2002c; Bell and 
Comeau 2003a). The intrinsic redshifts predicted 
by eqn 2 are listed in Table 1, for values less than 
z = 4.5. 

3. Results for the SDSS High-Redshift 
Data 

In models involving the ejection of QSOs from 
galaxies, in addition to the intrinsic component, 
Dopplcr components due to randomly directed 
ejection velocities and the Hubble flow will be 
present in the measured redshifts. An estimate of 
the size of the ejection velocities expected in such 
cases was obtained in the study of the 15 QSOs 
near NGC 1068 (Bell 2002a,c) where ejection ve- 
locities up to ~ 25, 000 km s _1 were found. When 
a cosmological component is added to this, even 
peaks separated by redshifts of z — 0.62 can easily 
be smeared out (Bell 2004). 

Equation 1 predicts that the spacing of the in- 
trinsic redshifts increases as the intrinsic redshifts 
increase. Thus the separation between the pre- 
dicted lines will be greater at high redshifts and 
are more likely to be free of the confusion that 
can occur at lower redshifts. Also, this redshift 
range above z = 2 has not previously contained 
a sufficiently high number of redshifts to allow a 
meaningful study. It now offers the opportunity 
to examine a completely different redshift range 
as well as a vastly more complete and indepen- 
dent sample, keeping in mind that there still may 
be selection effects in this portion of the data. It 
is important to realize as well, however, that selec- 
tion effects that miss sources can only affect those 
regions of the distribution where sources actually 
exist. Therefore it will not be easy to sort out 
exactly how much of a low-density valley in the 
redshift distribution is real and how much is due 
to the selection effect. 

For the quasars listed in the DR1 catalog, the 
redshift distribution between z = 2.3 and z = 5 is 
plotted in the lower histogram in Fig. 1. There 
is clearly a drop in the number of quasars at z = 
2.75. 

In this same redshift range the distribution of 
approximately 5000 quasar redshifts in the sample 
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Table 1 

Intrinsic Redshifts Predicted by Eqn 2 for z iQ < 4.5. 



(n) 


z zQ [N = l,n] 


z iQ [N = 2,n] 


z i(3 [7V = 3,n] 


ZiQ[N - 4, n] 


z«Q[iV = 5, n] 


z,q[N = 6,n] 


z iQ [JV = 7,n] 





0.620 


1.240 


1.860 


2.48 


3.10 


3.72 


4.340 
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Fig. 1. (upper histogram)Distribution of approximately 

5000 SDSS quasars between z — 2.3 and z — 5. (lower his- 
togram) Distribution of quasar redshifts in the DR1 catalog. 
Solid bars near the bottom indicate redshifts predicted by eqn 
1. Vertical dashed and dotted lines indicate redshifts predicted 
by eqn 2. 



containing 44,200 redshifts considered here (Oguri 
et al. 2003) is given by the upper histogram in Fig 
1. Most of the redshifts in this distribution have 
been obtained after several iterations of the target 
selection algorithm, and this process has clearly 
removed the narrow dip at z = 2.75. However, 
prominent peaks and valleys still remain. In fact, 
in some cases, the peak-to- valley ratio may be even 
higher than in the DM data. 

The locations of the preferred values predicted 
by eqn 1 are indicated in Fig 1 by short, bold ver- 
tical bars near the lower axis. There is clearly no 
agreement between these redshifts and the peaks 
in the distribution. This model is therefore not 
supported by these high-redshift data. 

The vertical dashed and dotted lines in Fig. 1 
indicate the locations of the intrinsic components 
predicted in this redshift range by eqn 2. Not only 
is there excellent agreement between these values 
and the peaks in the redshift distribution, there 
are also three clear valleys between the peaks. 
Valleys are expected here because no redshifts 
are predicted. They are not expected to drop to 
zero, however, because the preferred values will 
be smeared out by Doppler components. Not only 
are the peaks periodic in z = 0.62, they are pure 
harmonics above z = 0. While possible causes for 
the valleys at z = 2.7 and z = 3.5 have been sug- 
gested, no reason to suspect a valley at z = 4.03 
has been suggested. It is concluded here that eqn 
2 appears to explain the redshift distribution in 
Fig 1 better than the other arguments that have 
so far been proffered. Furthermore, although it is 
possible to suggest explanations for the valleys at 
z = 2.7 and 3.5, there is no way to prove that the 
reasons claimed are their true cause. If there are 
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actually no redshifts in the valleys, the selection 
effects raised above cannot affect the result. 

In Fig 1 it needs to be kept in mind from the 
magnitudes involved, that if these high redshift ob- 
jects have been ejected from nearby low-redshift 
galaxies, and are therefore still associated with 
them, they must be intrinsically several magni- 
tudes fainter than normal galaxies. Because they 
arc therefore too faint to be detected at large 
distances, their cosmological, or distance redshift 
components will be quite small relative to their 
intrinsic redshifts. In this model the distribution 
in Fig 1 is therefore essentially an intrinsic red- 
shift one. If, as is assumed in this model, quasars 
are ejected from active galaxies, and these galax- 
ies are distributed uniformly in space, the same 
should be true for quasars. Distant ones will not 
be detectable, however, with current sensitivities. 
On the other hand, those with lower intrinsic red- 
shifts, which have been found to be more lumi- 
nous than those with high intrinsic redshifts (Bell 
2002a, c), will be detectable to greater cosmologi- 
cal distances. 

4. The Full Sample of 44,200 Redshifts 

The entire sample containing 44,200 redshifts 
is plotted in Fig 2. As in Fig 1, from eqn 2, the 
vertical dashed lines indicate the highest redshift 
level in each TV-group and the dotted lines indicate 
the lower levels inside each group. Where the re- 
gion above z = 2.3 contains only two lines in each 
group, as can be seen in Fig. 2 and in Table 1, be- 
low z = 2 the density of predicted intrinsic redshift 
lines inside each group increases significantly. 

Immediately evident in Fig 2 is the fact that 
the density of SDSS redshifts in the lower half of 
the sample is much higher than in the upper half, 
although part of this difference can be attributed 
to the low-redshift sample with i-mag cut-off at i 
= 19.1. Similarly, there is a much higher density 
of predicted lines in the lower region, especially 
for the model described by eqn 2. Although, for 
eqn 2, the higher line density below z = 2 intro- 
duces some confusion, there is still evidence for 
peaks that correspond to the first three A-groups 
located at, or near, z = 0.62, 1.24, and 1.86, de- 
pending on the line distribution inside each group. 
Only at z = 2.48, where the source density is ris- 
ing steeply, is there no obvious peak, although the 
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Fig. 2. Distribution of 44,200 SDSS quasar redshifts. Solid 

bars near the bottom indicate redshifts predicted by eqn 1. 
Vertical dashed and dotted lines indicate redshifts predicted 
by eqn 2. 
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Fig. O. Correlation coefficient r between the distribution 

of 44,200 SDSS quasar redshifts and the distribution of lines 
predicted using eqn 2. The horizontal dashed line indicates the 
level above which r is considered significant. 
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DR1 distribution in Fig 1 does show one. The lack 
of a clear peak at z = 2.48 is not surprising since 
the survey portion with i-mag cut-off at i — 19.1 
found mainly sources below z ~2.2. This has con- 
tributed to the sharp drop in source density above 
z = 2.2, and such an abrupt change in level can 
easily mask a weak peak. 

In Fig 3 the correlation coefficient r, between 
the actual redshift distribution and the distribu- 
tion of rcdshifts predicted using eqn 2 is plotted 
versus shift. It shows a correlation coefficient of 
0.42 for zero shift. For a distribution containing 
45 samples as here, r greater than 0.288 is consid- 
ered significant. However, more importantly, as 
the two distributions arc shifted relative to each 
other, r first decreases and then increases again 
at ±6 bin shifts (one bin = 0.1 in redshift) corre- 
sponding to z — ±0.6, clearly indicating the pres- 
ence of peaks in the distribution that are periodic 
in z = 0.62. That the highest value of r is obtained 
for a shift of -0.6 is not surprising since it corre- 
sponds to the situation where the highest intrinsic 
line density (z ~ 1.2) is aligned with the highest 
peak in the quasar redshift distribution (z ~ 1.8). 
This further confirms that the periodic peaks in 
the distribution below z = 2 are also contributing 
to the correlation coefficient. The fact that r re- 
mains high over many shifts is due to the broad 
density correlation between the upper and lower 
halves of the data. 

There was no significant correlation found be- 
tween the SDSS redshift distribution and the dis- 
tribution of intrinsic redshifts predicted by eqn 1. 
This is not surprising since, in contrast to what 
was found for eqn 2, the intrinsic lines predicted by 
eqn 1 coincide with only one peak in the SDSS dis- 
tribution (z = 0.6), and there is little correlation 
between the source and line densities in the upper 
and lower halves of the distributions. It there- 
fore must be concluded that the model described 
by eqn 1 is not confirmed by the SDSS redshift 
data. Furthermore, while it has previously been 
demonstrated that the intrinsic redshifts predicted 
by eqn 2 also show a significant correlation with 
both the redshift distribution of the 574 strong 
quasars found in early surveys (Bell 2002c; Bell 
and Comeau 2003b), and the redshift distribution 
of the QSOs located around NGC 6212 (Bell and 
Comeau 2003b), there was also no significant cor- 
relation found in these cases using eqn 1. 



5. Conclusions 

In this examination of 44,200 SDSS quasar red- 
shifts it has been found that there are at least 5 
peaks in the redshift distribution that correspond 
to the pure harmonics of z = 0.62 predicted by eqn 
2. Although valid arguments have been put for- 
ward to explain the redshift valleys at z = 2.7 and 
3.5, there is no proof yet that they have correctly 
explained the causes of these low source counts, 
since eqn 2 predicts low densities at both these 
redshifts, reducing the effect of any selection ef- 
fect. Eqn 2 also predicts the low source count 
seen near z = 4.03 where no selection effect has 
yet been proposed. 

It is concluded here that the intrinsic redshift 
model given by eqn 2 can better explain the struc- 
ture present in the distribution of these 44,200 
SDSS quasar redshifts than can the selection ef- 
fects suggested by others. SDSS investigators are 
urged to keep this in mind if they hope to insure 
that real density fluctuations are not inadvertently 
removed. If we are to obtain a truly reliable sta- 
tistical sample, adjustments that affect the SDSS 
redshift distribution must not involve the assump- 
tion that quasar redshifts are purely cosmological. 
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